Quantum dynamics as a physical resource 
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I. INTRODUCTION 

The quantification and comparison of different types 
of physical resources lies at the heart of much of modern 
science. A good example is the physical resource energy, 
whose quantification enabled the development of ther- 
modynamics. More recently, motivated by applications 
to quantum information processing, there have been at- 
tempts to develop a quantitative theory of quantum en- 
tanglement [1]. This theory, still in its nascent stages, 
has been applied to gain insight into questions about the 
capacity of a noisy channel for information [2] , quantum 
teleportation with a noisy entangled resource [3] , and dis- 
tributed quantum computation [4]. 

Structurally, quantum mechanics has two parts, one 
part concerned with quantum states, the other with quan- 
tum dynamics. A general quantum dynamical process is 
described by a quantum operation (reviewed in [5]); such 
processes include unitary evolution, quantum measure- 
ment, dissipation, and decoherence. We believe quantum 
operations are a useful physical resource on an equal and 
logically independent footing to quantum states. 

The first step in studying a physical resource is to 
quantify it. Therefore, the purpose of our paper is to 
develop a theory quantifying the strength of quantum 
dynamical operations. Our motivations are axiomatic 
and operational questions concerning quantum dynam- 
ics. Our goal is to find strength measures capturing some 
of the structure in the complicated space of quantum 
operations, to gain insight into quantum dynamics and 
complex quantum systems [6, 7]. Although some of the 
measures we propose for operations are based on state 
entanglement measures, we expect the study of dynam- 
ics to provide different, complementary insights to those 
gained from the study of states. 

What questions will good strength measures allow us 
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to analyze? We foresee applications to the analysis of 
quantum computational complexity, distributed quan- 
tum computation, quantum communication, and quan- 
tum cryptography. As a simple example, consider the 
question of how many controlled-NOT (cnot) gates are 
required to implement a swap gate on two qubits, when 
assisted by arbitrary local unitaries. Suppose we have 
a measure K(U), quantifying the strength of a unitary 
U. Suppose further that K(U) satisfies (a) K{UV) < 
K(U) + K(V); and (b) K(U) = for local unitaries U. 
It is easy to see that the number of CNOT gates needed 
to do the swap gate is at least K (swAP)/if (cnot). 

More generally, the central problem of quantum com- 
putational complexity is to determine the minimum num- 
ber of one- and two-qubit gates necessary to implement 
a desired n-qubit unitary operation U . For example, U 
might encode the solution to a problem such as the travel- 
ing salesman problem. Suppose we have a strength mea- 
sure satisfying properties (a) and (b), above, as well as 
(c) K(U <8> I) = K(U). The number of gates needed to 
compute U is again bounded below by K (U)/K(cnot). 
Such a bound might help in determining the relation- 
ships between various quantum and classical complexity 
classes. We will return to this application several times. 

Another motivation to study quantum dynamics as a 
resource is recent work on universality in quantum com- 
putation. The class of interactions capable of performing 
universal quantum computation has been shown to be the 
class of bipartite entangling dynamics; any Hamiltonian 
which can create entanglement between any pair of qudits 
is universal, when assisted by arbitrary one-qudit uni- 
taries (see [8-13] and references therein, see also [14, 15] 
for related work). It has also been shown that any en- 
tangling two-qudit unitary, together with arbitrary one- 
qudit unitaries, is universal ([16], see [17] for a simple, 
constructive proof in the qubit case). 

These results show that there is a qualitative difference 
between entangling and non-entangling dynamics. Fur- 
thermore, they show all two-qudit entangling dynamics 
are qualitatively equivalent, as any one can simulate any 
other, provided local unitaries are available. By analogy 
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with the study of state entanglement, this suggests quan- 
tifying entangling dynamics. We now review prior work 
on this idea, organizing our discussion around three mo- 
tivating themes: the communication cost to implement 
an operation; the entangling ability of an operation; and 
the ability of an operation to communicate bits. 

The communication requirements for implementing a 
general bipartite unitary U were studied in Ch. 6 of [4] , 
where a general lower bound on the number of qubits of 
communication needed to implement U was proved, de- 
pending only on the operator Schmidt decomposition of 
U (see Sec. II for discussion of this decomposition). Eisert 
et al. [18] and Collins, Linden, and Popescu [19] studied 
the classical communication and entanglement required 
to implement some specific few-qubit quantum gates. 
Chcflcs, Gilson, and Barnett [20] studied the amount of 
communication and entanglement required to perform an 
arbitrary gate in a network of qubits. 

The capacity of a quantum operation to gener- 
ate entanglement seems to have first been studied by 
Makhlin [21], who found three invariants characterizing 
the non-local properties of two-qubit unitaries. Makhlin 
used the invariants to obtain results about entanglement 
generation, with a view towards applying them to the 
complexity of implementing gates. Zanardi, Zalka, and 
Faoro [22], Zanardi [23], and Wang and Zanardi [24], all 
obtained results about the average entanglement gener- 
ated by a unitary. Cirac et al. [25] studied the ability 
of an operation to produce entanglement by mapping 
the operation onto a corresponding state, and studying 
the properties of that state. Kraus and Cirac [26] stud- 
ied the maximum entanglement which can be created by 
a unitary operator acting on two initially unentangled 
qubits. They found an explicit formula for the maxi- 
mum entanglement that can be generated without ancil- 
las, and showed that this amount can be exceeded with 
the use of ancillas. Leifcr, Henderson, and Linden [27] 
used similar reasoning to obtain an explicit formula for 
the entanglement generated without ancillas, but allow- 
ing initial entanglement. They also obtained numerical 
results demonstrating that the addition of ancillas can 
increase the maximum entanglement generated. In a dif- 
ferent context, Scarani et al [28] related the entangling 
power of a unitary operation to the problem of thermal- 
ization of a quantum system. 

A related approach is to quantify the entangling abili- 
ties of Hamiltonians rather than unitaries. Diir et al. [10] 
considered the rate at which a Hamiltonian creates en- 
tanglement, and found techniques to optimize this rate. 
More recently, Vidal, Hammerer, and Cirac [29] (see 
also [30]) analytically characterized the minimum time 
required to simulate one Hamiltonian with another, and 
found the minimum time required to simulate a desired 
unitary with a Hamiltonian. This allowed them to de- 
fine a partial order on unitaries, according to which one 
unitary U is more non-local than another unitary V if 
and only if, for any Hamiltonian, the minimum time re- 
quired to simulate U is longer than the minimum time to 



simulate V. They also obtained results on the optimal 
choices of non-local interactions for transmitting classical 
bits between two parties. Childs et al. [31] found an ex- 
plicit formula for the maximum entanglement created by 
a class of two-qubit Hamiltonians, including the Ising in- 
teraction and the anisotropic Heisenberg interaction, for 
which this maximum is achieved without ancillas. 

The ability of a quantum operation to communicate 
classical information was studied by Beckman et al. [32] , 
who obtained simple necessary and sufficient conditions 
for information transmission to be possible. Bennett et 
al. [33] and Berry and Sanders [34] studied the capacity 
of a bipartite operation to communicate information, and 
related this capacity to the ability of the operation to 
generate entanglement. 

Our paper draws on all these perspectives, but differs 
in an important way. Rather than focusing on the abil- 
ity of a quantum dynamical operation to generate some 
static resource, such as entanglement or shared classical 
bits, we believe it is possible to quantify quantum dy- 
namical operations as a physical resource in their own 
right. That is, we do not need to make reference to the 
ability of the operation to generate some other resource. 

How can one develop a theory of dynamic strength 
without relying on familiar state-based resources? The 
approach we take is to identify plausible axioms and 
properties a good measure of strength should satisfy, and 
develop measures satisfying those properties. 

The paper is structured as follows. Sec. II opens 
by introducing two concrete examples of strength mea- 
sures for unitary operations, the Hartley and Schmidt 
strengths. Sec. Ill considers operational questions moti- 
vating strength measures, and uses these questions to mo- 
tivate some abstract axioms for such measures. Sec. IV 
briefly summarizes a useful canonical decomposition for 
two-qubit unitary operators. Secs.V and VI explore a 
variety of specific definitions for dynamic strength mea- 
sures. Our general philosophy is to explore a wide vari- 
ety of measures, and then to concentrate on those which 
appear most likely to yield useful practical answers to in- 
teresting operational questions. Sec. VII concludes with 
a summary and a table of results. 



II. INVITATION: THE HARTLEY AND 
SCHMIDT STRENGTHS 

In this section we introduce two strength measures, 
the Hartley strength and the Schmidt strength. These 
measures are introduced both because of their intrin- 
sic interest, and also because the Hartley strength will 
be used as a simple, concrete example illustrating the 
more abstract, axiomatic approach to dynamic strength. 
The Hartley and Schmidt strengths are based on a gen- 
eralization of the Schmidt decomposition to operators, 
which we call the operator-Schmidt decomposition [4]. 
To explain the operator-Schmidt decomposition we in- 
troduce the Hilbcrt- Schmidt inner product on the space 
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of d x d operators, defined by (Q,P) = tr(Q^P), for 
any operators Q and P. Using this inner product, we 
define an orthonormal operator basis to be a set {Qj} 
which satisfies the condition: (Qj, Qk) = tr(Q^Qk) = 5jk 
where Sjk = 1 if j = k and otherwise. For exam- 
ple, a complete orthonormal basis for the space of one- 
qubit operators is the set of normalized Pauli matrices, 
{/, X, Y, Z} = {//A X/y/2, Y/y/2, Z/V2}. 

An operator Q acting on systems A and B may be 
written in the operator-Schmidt decomposition [4]: 



(2-1) 



where si > and Ai and Bi are orthonormal operator 
bases for A and B, respectively. To prove the operator- 
Schmidt decomposition, expand Q = X^fc MjkCj <S> D^, 
where Cj and Dk are fixed orthonormal operator bases 
for A and B, respectively, and Mjk are coefficients. The 
singular value decomposition states that the matrix M, 
whose (j, fc)th entry is Mjk, may be written M = UsV, 
where U and V are unitary, and s is diagonal with non- 
negative entries. We thus obtain 



jkl 



D k , 



(2.2) 



where si is the Ith diagonal entry of s. Defining Ai = 
y\- UjjCj and Bi = ^frVikDk, which are easily shown 
to be orthonormal operator bases for A and B, we obtain 
the operator-Schmidt decomposition Eq. (2.1). 

Nielsen [4] defines the Schmidt number of an operator, 
Sch(Q), to be the number of non-zero coefficients in the 
operator-Schmidt decomposition for Q [68]. 

A simple example is the CNOT gate which has operator- 
Schmidt decomposition 



CNOT 



= V2\0) 



y/2\l)(l\®X 



(2.3) 



and hence has Schmidt coefficients {\/2, V%}, and 
Sch(CNOT) = 2. The swap gate for qubits has operator- 
Schmidt decomposition 



SWAP =I®I+X®X+Y®Y+Z®Z (2.4) 

and hence Sch(swAp) = 4. A less familiar example is the 
gate 



U p = (yjl-pl <g> J + iy/pX (g>X^jx 



l-pl ® J + i s /pZ ® Z 



which has operator-Schmidt decomposition 



(2.5) 



U p = 2 (1 - p) I ® I + 2pY O Y+ 

2y/p(l-p) [(e 4 ^ 4 !) ® X + (e iw / 4 Z^j ® Z 

(2.6) 

and thus has Schmidt number 1 when p = or 1 , and 4 
otherwise. 



A more complicated example is provided by the quan- 
tum Fourier transform, whose unitary action on I qubits 
is defined by the action on computational basis states [35] 



^ 2-1 



\t), 



(2.7) 



where we number the basis states from |0) through 
1 2 l — 1). A useful alternate formula for the quantum 
Fourier transform may be obtained by working in a bi- 
nary representation, s = s\ ■ ■ ■ si, whence [69] 
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l/.,>®l/ a ,_ iai > 



I/ S1 - 



(2.8) 



where the one-qubit state \f t ) is defined for an arbitrary 
bit string* = ti • • -t k by [//) = [|0)+exp(27ri0.i)|l)]/A 

and O.t is the binary fraction ti/2 + t%/4 -\ + t k /2 k . 

Suppose now that system A consists of m qubits and 
system B consists of n qubits, and U is the quantum 
Fourier transform on to + n qubits. From Eq. (2.8), 

U\xi,...,x m ,yi,...,y n ) = \ fy n )®---®\f Xl - yn )- (2.9) 

Suppose to < n. To determine the Schmidt decom- 
position of the quantum Fourier transform it is conve- 
nient to introduce the notation y' = yi ■ ■ ■ y n -m an d 
y" = j/n-m+i • • 'Vni so the string y can be formed by 
concatenating the strings y' and y" . It follows from the 
previous equation that 



U 



xy" 



B 



xy" i 



(2.10) 



where x ranges over m-bit strings x\ ■ ■ ■ x m , and we define 

A xy ii = \ fy n )\fy n - 1 y v ) ' ' ' \ fy n - m+ i---y n ) { x \ 



B. 



xy" 



a 



xy'y" 



7 j Cxy'y" 

y' 

\fyn-m---Vn 



(2.11) 



\fx 1 - Vn ){y\- 



A calculation shows that the A xy n are orthonormal 
operators, and the B xy " are an orthogonal set, with 
(B xy ii , B xy ii) = 2™~"\ Thus the Schmidt decomposition 
for the quantum Fourier transform is 



U 



xy" 



B 



xy" 



(2.12) 



Thus, when to < n the quantum Fourier transform has 
Schmidt number 2 2 m , and all nonzero Schmidt coeffi- 
cients are equal to v / 2™~™. Note that the Schmidt de- 
composition of the quantum Fourier transform was al- 
ready obtained in [4] when to = n; we have not yet suc- 
ceeded in determining the Schmidt decomposition of the 
quantum Fourier transform when to > n, but conjecture 
that it has Schmidt number 2 2 ™[70] 

The Hartley strength [71] of an operator K^ a , T (Q) is 
defined by 



K Sai (Q) = log[Sch(Q)]. 



(2.13) 
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(The logarithm is taken to base 2 throughout this paper.) 
Returning to our examples, the CNOT gate has Hartley 
strength log 2 = 1, the SWAP gate has strength 2, and 
Up has strength for p = or 1, and strength log 4 = 2 
otherwise. The quantum Fourier transform has Hartley 
strength 2m, provided m < n. 

The Schmidt strength is motivated by a simple obser- 
vation about unitary operators U acting on systems A 
and B of respective dimensions dj± and d®. For such an 
operator, the relation tr(WU) = d^d® implies that the 
Schmidt coefficients s ; satisfy J2i s f = ^A^b- Therefore, 
the numbers sf / (d^ds) form a probability distribution. 
A natural measure of the non-local content of U is thus 
the Schmidt strength, defined to be the Shannon entropy 
H(-) of the distribution sf/(<Ucfe) [72], 

X Sch ([/)EEff({^}). (2.14) 

More generally, for an arbitrary bipartite operator Q we 
define the Schmidt strength by 

K Sch (Q) = H ({to(QTQj})> (2-15) 

where {sf /tr(Q^Q)} are the squared Schmidt coefficients 
of Q, normalized to form a probability distribution. Note 

that #Sch(CNOT) = 1, if Sch (SWApj = 2, K Sch (U p ) = 
H[(l-p) 2 ,p 2 ,p(l-p),p(l-p)], and K Sch = 2m for 
the quantum Fourier transform, when m < n. 

III. CONCEPTUAL FRAMEWORK 

In this section, we explore two approaches to the defi- 
nition of strength measures. In the operational approach, 
discussed in Sec. Ill A, we define several measures of 
strength based on the ability of an operation to perform 
various tasks. These measures thus quantify a dynami- 
cal resource required by each task. The second approach, 
the axiomatic approach, is explored in Sec. Ill B, where 
we identify a list of three axioms and nine useful proper- 
ties for a strength measure. These two approaches may 
appear to be independent, but there is actually substan- 
tial interplay. In particular, many of the properties in 
Sec. Ill B are motivated by consideration of the opera- 
tional measures of strength of Sec. Ill A. 

A. Operational approach 

Quantum dynamics are clearly an essential component 
in quantum information processing tasks. However, it is 
difficult to identify which properties of quantum dynam- 
ics are the most essential, because different properties 
are required for different tasks. This variety is reflected 
in this section by the fact that different operational ques- 
tions give rise to different notions of strength. 

The reader should note that the main point of this sec- 
tion is not to prove results about the measures we define. 



Rather, it is to provide definitions of some strength mea- 
sures, and a discussion of their operational motivation. 
After we have enumerated the properties we would like 
these measures to satisfy, we will take up the problem 
of determining the properties of these measures, and the 
relationships between them. 

1. Entanglement generation and communication capacity 

In this section, we consider two related questions on 
the ability of a quantum operation to create entangle- 
ment and to communicate information. We also review 
some of the recent work on these subjects. 

How much entanglement can be generated by a 
quantum operation? 

How much entanglement a single application of a uni- 
tary U can generate depends crucially on the initial states 
U may act on. We must also specify whether we are in- 
terested in the maximum, minimum, or average entangle- 
ment generated. We focus primarily on maximizations. 

We define two measures for the entangling strength of 
a unitary U. (See Sec. V A3 for some generalizations to 
quantum operations.) 

The first strength measure quantifies the maximum en- 
tanglement which a unitary U can create between two 
systems A and B with the use of arbitrary ancillas, but 
without prior entanglement: 

K E (U)=mHxE{U\a)\fl) (3.1) 

\c),\0) 

where |a) ranges over all (possibly entangled) states of 
system A plus an ancilla TZa, and \(3) ranges over states of 
system B plus an ancilla 1Z&, and E is the usual measure 
of bipartite pure state entanglement, the von Neumann 
entropy of the reduced density matrix [73]. Note that 
the ancillas may be chosen with dimensions equal to the 
dimensions of A and B respectively, since the Schmidt 
number of \a) with respect to the A : 1Zj± division is at 
most d_A, and similarly for \j3). It follows that Ke is truly 
a maximum, and not a supremum. 

Kraus and Cirac [26] calculated Ke{U) for some spe- 
cial two-qubit unitaries, while Leifer, Henderson and Lin- 
den [27] obtained numerical evidence that removing the 
ancillas decreases the maximum entanglement for certain 
unitaries. 

The second measure allows the possibility of prior en- 
tanglement as well as ancillas. Kae(U) is the magnitude 
of the maximal change in entanglement caused by U: 

K AE (U) = sup - E(\m , (3.2) 

IV-) 

where \ip) ranges over all states of ATZ^ and BIZb [74]. 

Clearly, Ke{U) < Kae(U) for all U. Later, we will see 
that there exist unitaries U for which Ke{U) ^ K\e{U), 
demonstrating that these two measures capture different 
notions of a unitary 's ability to generate entanglement. 
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An alternative approach to quantifying entanglement 
generation has been explored by Zanardi [23] , and Wang 
and Zanardi [24]. Zanardi [23] defines a measure of en- 
tanglement, L(U), for a unitary operator U onarf^xdg 
system by the linear entropy, L(U) = 1 — s\jd\d 2 B , 
where si are the Schmidt coefficients of U. Provided 
d-A = d-B = d, it can be shown that [23, 24], 



/ 



dad/3 L[U(a®0)] 
d 2 



(3.3) 



[L(U) + L([/SWAP) - L(SWAP)] 



where da and dj3 are the uniform, normalized, Haar mea- 
sures on the first and second qudits, respectively, the 
function L on the left is the measure of state entangle- 
ment based on the linear entropy of the squared Schmidt 
coefficients of the state, while the function L on the right 
is the operator entanglement defined by Zanardi. This 
equation nicely connects the Schmidt coefficients and the 
average entanglement generated by U. 

In a similar vein, Wang and Zanardi [24] define a notion 
of concurrence for unitary operators with Schmidt num- 
ber 2. For a system AB of dimension d^xdjs, they define 
C(U) = 2s\S2/{dj(dB), where s\ and S2 are the Schmidt 
coefficients of U . This definition extends the notion of 
concurrence for qubits introduced by Hill and Woot- 
ters [36]. Simple algebra and the fact that Y^i s f = dji.de 
implies that C 2 (U) = 2L(U), where L(U) is the measure 
of operator entanglement introduced by Zanardi [23] . 
How useful is a quantum operation for communi- 
cation? 

An interesting question is to determine the relation- 
ship between the entanglement generated by a channel 
and its capacity to transmit classical information between 
two systems. Recently, Bennett et al. [33] and Berry and 
Sanders [34] have examined the relationship between the 
entangling capacity of a two-qubit unitary and its abil- 
ity to transmit information. In particular, Bennett et al. 
considered the maximum entanglement that can be gen- 
erated from any (possibly entangled and mixed) state 
with t uses of the unitary gate U. They argued that the 
maximum entanglement generated with t uses of U is just 
t times the maximum entanglement generated with one 
use of U, and that Kae is an upper bound on the av- 
erage number of bits which can be reliably transmitted 
between A and B. 



2. Quantum computational complexity 

In this section we consider a different motivation for 
the study of quantum dynamics as a resource. Rather 
than considering an operation's explicitly non-local 
properties (such as its ability to create entanglement), 
we ask what characterizes the difficulty of performing a 
quantum computation. 



A reasonable measure of the complexity of implement- 
ing a unitary U with a gate set U is simply the min- 
imum number of gates from U in a circuit which im- 
plements U. For example, suppose we only have the 
ability to implement the CNOT gate on two qubits, with 
either acting as the control, and we wish to simulate 
the SWAP gate. In this case we have the gate set 
U = {CNOT12, CNOT 2 i} where the first subscript refers 
to the control qubit and the second the target. Since 
SWAP = CNOT12CNOT21CNOT12 (and the SWAP gate can- 
not be implemented with only two CNOT gates) , the com- 
plexity of the swap gate relative to U is 3. 

To generalize this idea, we define K com : 

K com (U\V) = min . x (Wj) \u = J] W jt Wj G U } , 

(3.4) 

where the cost function x(Wj) IS an y non- negative func- 
tion that quantifies the difficulty associated with imple- 
menting Wj. 

The circuit complexity measure has the property that, 
for any two unitary operators U and V, 

K com (UV\V) < K com (U\V) + K com (V\V), (3.5) 

since one circuit implementing UV is the concatenation 
of the minimal circuits implementing V and U separately. 
We refer to this property as the chaining property. 

In general, K com is prohibitively difficult to calculate 
since it is very hard to prove that a given circuit for 
U is minimal. However, it is possible to find lower 
bounds on K com as follows. Expanding upon the ex- 
ample given in the introduction, suppose U is a two- 
qudit unitary, and one is given the ability to perform 
a set of two-qudit gates U = {U\, . . . ,U m }, and lo- 
cal unitary operations. What is the minimum number 
of two-qudit gates required to implement U? Suppose 
U = (A <g> BQ)U h {Ax <g> B x ) ■■■U h (A k <g> B k ), where 
Aj (g) Bj denotes a local unitary, and Ui j G U. Let K 
be any measure satisfying K(UV) < K{U) + K{V) and 
K(A <g> B) = for any local unitary A® B. Then 

K{U) = K[(A ® B )U h (Ax <g> ■ ■ ■ Ui„ (A k ® B k )] 
<K(U h ) + --- + K(U lk ) 



< kK„ 



(3.6) 



where K max is the maximum value of KiU^). We have 
deduced a useful bound on the number of gates, 



k > 



Kim 

K\n ax 



(3.7) 



This captures the intuitively appealing notion that the 
number of gates required to implement U is at least 
equal to the total strength of U, divided by the maxi- 
mum strength of any of the implementing gates. Indeed, 
if we take the cost of a local unitary to be and the 
cost of a two-qudit gate to be 1, the argument implies 
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that K com (U\V) > K(U)/K maK . Although this argu- 
ment holds only for two-qudit unitaries, U, we will ex- 
tend it to n-qudit unitaries after the discussion of stabil- 
ity properties in the next section. 



B. Axiomatic approach 

One approach to quantifying entanglement is to con- 
sider axioms which an entanglement measure "ought" 
to satisfy, and to explore the consequences of those ax- 
ioms [1, 37-39]. While this approach has occasionally 
been criticized [40] , it has certainly proven fruitful. 

Here we explore an analogous axiomatic approach to 
the study of strength measures for quantum dynamical 
operations. We propose a number of axioms that such 
measures might be expected to satisfy and investigate 
some implications of these axioms [75]. 

The structure we adopt is to first describe (in 
Sec. Ill B 1) the fundamental axioms that we expect any 
strength measure should satisfy. We then describe some 
other useful properties a strength measure may satisfy in 
Sec. Ill B 2. Finally Sec. Ill B 3 illustrates the axiomatic 
framework by applying it to the analysis of the commu- 
nication cost of distributed quantum computation. 



1. Fundamental properties 

We denote our strength measure by K(£), where £ is 
a trace-preserving quantum operation acting on a set of 
n systems, A\, . . . ,A n , of dimensions d\, . . . ,d n . We will 
frequently be interested in the case where £ is a unitary 
quantum operation £{p) = UpU^ for some unitary U. In 
this case, we write K(U) to denote the dynamic strength 
of U. We will also use the convention that the symbol 
for a unitary such as U may either mean the unitary op- 
erator U , or the corresponding quantum operation, that 
is U(p) = UpU^. This abuse of notation will only be 
employed when its meaning is clear from context. 

As each axiom is introduced we illustrate it by exam- 
ining whether the Hartley strength satisfies the axiom. 
Note that Kn av (U) is defined for a unitary operator U 
acting on two systems labeled A and B of dimension dj± 
and dg, respectively. 

Axiom 1 (Non-negativity) K{£) > for all quantum 
operations £. 

This is more a convention than an axiom, which we intro- 
duce as a convenience to simplify many of the properties 
below. The Hartley strength satisfies this axiom. 

Axiom 2 (Locality) K(U) > with equality if and 
only if U can be written as a product of local unitary 
operations. 

The Hartley strength Kn ar (U) satisfies locality. 



The axiom of locality captures the idea that dynamic 
strength measures the non-local content of a quantum 
gate. For example, in the bipartite case, it is possible to 
generate entanglement with a unitary U if and only if U 
cannot be written as a product of local unitary opera- 
tions. Similarly, it is possible to communicate classical 
information with a unitary if and only if it cannot be 
written as a product of local unitaries [32]. Summariz- 
ing, for any K satisfying locality, we have K(U) > if 
and only if U is capable of generating entanglement or, 
alternatively, of transmitting classical information. 

How should the axiom of locality be extended to non- 
unitary operations? For example, we might require that 
K(£) > if and only if £ cannot be implemented by local 
operations and classical communication. Or perhaps we 
might require that K{£) > if and only if £ generates 
quantum states with non-zero entanglement (according 
to some entanglement measure). Many other possibilities 
can be imagined which we will not enumerate. 

Axiom 3 (Local unitary invariance) Suppose 

A\ , . . . , A n and B\ , . . . , B n are local unitary operations 

on the respective systems A\ , . . . , A n . Then 

K [(Ax ® • • • ® A n ) o £ o (Bi o • • • o B n )] = K{£). 

(3.8) 

The Hartley strength satisfies local unitary invariance. 

The axiom of local unitary invariance requires that the 
strength of a quantum operation is not changed by local 
operations. Thus, it is in accord with the notion that the 
strength is a measure of an operation's non-local content. 

2. Other useful properties 

We have just introduced three axioms essential for any 
strength measure describing the non-local content of an 
operation. We now introduce several useful properties a 
strength measure may satisfy, beginning with two invari- 
ance properties. 

Property 1 (Exchange symmetry) Let £ be a quan- 
tum operation acting on a multipartite system whose sub- 
systems have the same Hilbert space. The swap opera- 
tion acting on any two of these components has the effect 
of interchanging their states. Then K has the exchange 
symmetry property if for all such SWAP operations, 

if(swAP o £ o swap) = K{£). (3.9) 

Property 2 (Time-reversal invariance) For all uni- 
taries U, K(W) = K{U). 

The Hartley strength satisfies both axioms. 

Property 3 (Continuity) For some metric D(-, •) on 
the space of quantum operations, \K{£) — K{T)\ < 
f(D{£,!F)), where f{-) is a continuous and monotone 
increasing function such that /(0) = 0. 
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The Hartley strength is not continuous with respect to 
standard metrics on the space of unitary operations: the 
presence of any non-locality in a unitary operation U is 
sufficient to cause a discontinuous jump in the Hartley 
strength from to 1 or more. 

A major use of the continuity property is in the anal- 
ysis of quantum computational complexity problems; see 
the discussion after the chaining property 

Property 4 (Chaining) Suppose £ and T are two 

quantum operations. Then K{£ o T) < K{£) + K{F). 

The main utility of chaining was anticipated in the 
introduction: it can give bounds on the number of gates 
required to perform a particular quantum operation. 

When combined with the continuity property the 
chaining property may also be used to prove bounds 
on the approximation of unitary operations. This is 
important in applications to computational complexity 
since it is usually sufficient to solve problems with a 
high probability of success. Suppose, for example, that 
U is a desired two-qudit unitary operation, and one is 
given the ability to perform a set of two-qudit gates 
U = {Ui, . . . , U m }, and local unitary operations. Let 
K be any measure satisfying continuity, for some choice 
of / and D, as above, as well as chaining and locality. Let 
Aj (g> Bj be local unitaries and Ui j G U. To obtain an ap- 
proximation V = (A (g> B )U h (A 1 (g> Bi) ■ ■ ■ Ui k (A k ®B k ) 
to U such that D{U, V) < e we need, by the continuity 
property K(V) > K{U) - /(e). But K(V) < kK max , 
where K maK is the maximum value of K (Ui), so the num- 
ber of gates satisfies 



k > 



k(u) m 



K„ 



(3.10) 



The Hartley strength satisfies the chaining property, 
but to prove it we need a related lemma. 

Lemma 1 Suppose U has operator- Schmidt decomposi- 
tion U — ^2jSjAj ® Bj. Suppose U can be written in 

some other form as a sum over products, U = ^ fe A k <8> 
Bk ■ The number of terms in this decomposition is at least 
as great as the number of terms in the operator- Schmidt 
decomposition. Thus, the operator- Schmidt decomposi- 
tion is a minimal decomposition for U , in the sense that 
it has the fewest product terms of any sum- over-products 
decomposition. 

Proof: A simple proof of the lemma is to note that: 
Aj = -ti B \(I ® B))U] =-J2 4fctr(Bjs fc ). (3.11) 

Thus each Aj can be written as a linear combination of 
the Ak- But the Aj are orthonormal, and thus linearly 
independent. It follows that the number of operators A k 
must be at least as great as the number of Aj, that is, at 
least as great as the Schmidt number of U. ■ 



With Lemma 1 in hand it is straightforward to prove 
the chaining property. Suppose U = s jAj ® Bj and 
V = J2k tkCk ® Dk are Schmidt decompositions for uni- 
tary operators U and V. Then we have 



UV = J2 Sjtk(AjC k ) ® {BjD k ). (3.12) 

jk 

The total number of terms in this sum-over-products 
decomposition of UV is Sch([/)Sch(F), and so by the 
lemma we must have Sch(UV) < Sch(V)Sch(F). Tak- 
ing logarithms of both sides of this inequality yields the 
chaining property for the Hartley strength. 

Until now we have only been concerned with strength 
measures defined for fixed quantum systems. Compare 
this with the situation for entanglement measures. It is 
often said that there is a unique [39, 41, 42] entanglement 
measure for bipartite pure states, namely, the von Neu- 
mann entropy of the reduced density matrix. Strictly 
speaking, this is not a single entanglement measure, since 
it can be applied to many different types of quantum sys- 
tems — pairs of qubits, a qubit and a qutrit, and so on. 
Rather, it is a family of entanglement measures, satisfy- 
ing certain consistency properties that make it sensible 
to refer to it as a single measure. 

Motivated by this, we describe two consistency prop- 
erties we expect of a family of strength measures. There 
are two different ways in which a family of strength mea- 
sures arises naturally. The first corresponds to appending 
additional systems while keeping the state-space dimen- 
sions of the existing systems constant. The second cor- 
responds to fixing the number of systems, and varying 
the state-space dimensions of the individual systems by 
adding local ancillas. 

For the statement of each of the following properties 
we imagine that there is a family of strength measures, 
each of which is denoted by the same letter K. When 
necessary, we add superscripts to make precise which sys- 
tems K is acting on. For example, K A ' B ' C {£) indicates 
the strength with respect to a division into three com- 
ponents, labeled A, B, and C, and K A ' BC (£) indicates 
the strength with respect to a division into two compo- 
nents A and BC. For notational simplicity, we state these 
properties for the case of three systems, with the gener- 
alization to more systems following similar lines. 

Property 5 (Stability under addition of systems) 

Suppose £ acts on systems A and B, and C is an addi- 
tional system. Then the family K is stable with respect 
to additional systems if 



K A:a {£) > K A {£ (8) J), 
where I denotes the identity operation on C. 



(3.13) 



Note that it does not make sense to speak of the Hartley 
strength as being stable or not stable in this sense, since 
it is only defined for two-component systems. 

The intuition motivating the inequality in the state- 
ment of stability is that the "two-party" non-locality 



8 



present in £ should not be less than the "three-party" 
non-locality in £ ®I. A stronger statement of the stabil- 
ity property would replace the inequality by an equality. 

The stability property is useful in the context of quan- 
tum computational complexity. We explained earlier how 
to derive lower bounds such as Eq. (3.7) and Eq. (3.10) 
on the number of gates needed to implement a two-qudit 
quantum operation. In the context of quantum compu- 
tational complexity, the most natural setting is that we 
wish to implement a family of n-qubit unitaries U (in- 
dexed by n) using a universal set of one- and two-qubit 
quantum gates. In such a setting, we are looking for the 
most efficient decomposition of each U into a product of 
two-qubit gates 

U = U jlkl U hk2 ---U jlkl , (3.14) 

where the subscripts denote the qubits on which each 
(possibly different) unitary gate acts. A bound on the 
minimum number of gates I may be deduced from the 
chaining and stability properties, using a similar anal- 
ysis to that given in connection with chaining alone, 
I > K(U) I 'iimax, where now K max is the maximum value 
of the strength of any two-qubit gate. Because of sta- 
bility, -Kmax is a constant, independent of n, so in order 
to prove interesting lower bounds on I, one only needs 
to analyze the asymptotic behavior of K(U) as a func- 
tion of n. If, for example, we could find a strength mea- 
sure satisfying both chaining and stability, and such that 
K(U) = 9(2") for some family of unitaries U, then it 
would follow that the family requires a number of gates 
exponential in n. If, in addition, K has suitable continu- 
ity properties, then it may be possible to prove that the 
family requires exponential time even if some reasonable 
probability of error is allowed. Needless to say, if this 
were true for a unitary encoding of, say, the solution to 
a problem such as the traveling salesman problem, this 
would be a very interesting result indeed. 

Our second notion of stability is that introducing local 
ancillas which arc then ignored should not change the 
strength of an operation. 

Property 6 (Stability with respect to ancillas) 

Suppose £ acts on systems A and B, and C is an 
additional system. Then the family K is stable with 
respect to local ancillas if 

K A -- B (£) = K A:BC (£(g>l). (3.15) 

The Hartley strength is clearly stable with respect to 
local ancillas. 

We now move on to additivity properties. 

Property 7 (Weak (sub)additivity) Suppose Ai, 
Ai, B\, and B2 are distinct systems such that A\ and 
A2 have the same state space, as do B\ and B2- Suppose 
£ is a quantum operation that can act on either A\B\ 
or A2B2- Then the family K is weakly subadditive if 

K MM '- B ^{£®£) <2K Ai:Bl (£). (3.16) 



K is weakly additive if the inequality can be replaced by 
an equality in the above expression. 

Property 8 (Strong (sub)additivity) Suppose A\, 
A2, B\, and B2 four distinct systems, and £ and T 
are quantum operations acting on A\B\ and A2B2, 
respectively. Then the family K is strongly subadditive 

if 

^■■ Aj:BlB! {£ ® T) < K Al ' Bl {£) + K A ^ [T). (3.17) 

K is strongly additive if the inequality can be replaced by 
an equality in the above expression. 

Note that strong subadditivity for a strength measure is 
not connected with the strong subadditivity property for 
quantum mechanical entropy [43] . 

The Hartley strength is strongly additive for unitary 
operations U and V, and thus possesses all four of these 
properties. To see this, suppose U and V are unitary 
operators with Schmidt decompositions U = ^ . SjAj ® 
Bj and V = J2 k tkC' k ® D k , where Aj, Bj, C k , and D k 
act on systems A\,B\, ^2, and B2, respectively. Then the 
Schmidt decomposition of U ® V with respect to A1B1 : 
A2B2 is 

U <g> V = s M A i ® C fc ) <8> {Bj ® Li fe ). (3.18) 

jk 

It follows that Sch(f7 <g> V) = Sch([/)Sch(V) and, taking 
logarithms, we see that the Hartley strength is strongly 
additive. 

Proposition 2 // the family K satisfies the chaining 
property and is stable with respect to local ancillas, then 
it is strongly subadditive. 

Proof: Applying simple algebra, the chaining property, 
and stability with respect to local ancillas in turn, we 
have: 

K(£ <g) T) = K[{£®1){1®T)] 

<K{£®1) + K{1®T) (3.19) 
<K{£)+K{F), 

which is the strong subadditivity property. ■ 
The converse is not true — we will see later that the 
Schmidt strength is strongly additive and stable with re- 
spect to local ancillas, but does not satisfy chaining. 

The final property addresses what happens when a 
quantum operation arises as a consequence of tracing out 
part of the action of a quantum operation acting on a 
larger system. For notational simplicity, we state this 
property for the special case of two systems, with the 
generalization to more systems following similar lines. 

Property 9 (Reduction) Suppose a quantum opera- 
tion £ on a composite system AB is obtained from a 
quantum operation on ABC as follows: 

£(pAB)=trc{F( P AB®o- c )}, (3.20) 
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for some fixed state o~c °.f system C. Then a family 
K of strength measures has the reduction property if 
K A:B (£) < K A:BC (F). 

The intuition behind the reduction property is that if it is 
possible to do T , then it is also possible to do £, without 
any extra dynamical resources being required. 

The reduction property is important both in the anal- 
ysis of distributed quantum computation (see below) and 
for the applications to quantum computational complex- 
ity suggested earlier in this paper. In the latter applica- 
tions we implicitly assumed that the implementation of 
some desired unitary could not be assisted by the in- 
troduction of ancilla qubits that are discarded at the 
end of the computation. However, there is evidence to 
suggest that ancilla may help in performing a unitary 
transformation quickly; for example, some of the con- 
structions in [44] were made more efficient by the use 
of ancilla. Suppose, however, that K has the reduction 
property, and that U can be implemented by perform- 
ing an operation Vona larger system. That is, suppose 
V|^)|s) = (U\tp))\s'), for all \ip), and for some fixed an- 
cilla states \s) and |s'). Then we have K(U) < K(V). 
If, in addition, it is possible to use K(-) to prove bounds 
on computational complexity, as described earlier, then 
it follows from the inequality K(U) < K(V) that any 
bound on the computational complexity of U must also 
apply to V, and thus our techniques can be applied even 
when working qubits are allowed. 

The reduction property makes restricted sense for the 
Hartley strength, which is defined only for unitary oper- 
ators. In particular, imagine, as above, that we have a 
unitary V acting on ABC such that V\ip)\s) = (U\ip))\s'), 
where is an arbitrary state of AB, U is a unitary act- 
ing on AB alone, and \s) and \s') are fixed states of C. 
To see that i^Har satisfies the reduction property, let us 
introduce orthonormal bases \j), \k), and |Z) for the sys- 
tems A, B, and C, respectively. Note that the invariance 
of i^Har with respect to unitaries on system C implies that 
it suffices to consider V such that y|-i/>)|0) = U\ip)\0), 
where |0) is the first element of the basis for C. Suppose 
we expand V as 

V= Yl V j , k , l >, jk i\f)(j\®\k'l')(kl\, (3.21) 

jklj' k'V 

where the comma in the subscript of V separates the 
row index from the column index. Since and 
\k'l')(kl\ arc orthonormal operator bases, it follows that 
the Schmidt coefficients of V are just the singular val- 
ues of the matrix V defined by Vjj\kk'U' = Vj'k'l',jkl- 
Thus, the Schmidt number of V is given by the number 
of non-zero singular values, or the rank, of the matrix . 
Similarly, we can expand U as 

U= U j'k<,jk\3 , )(j\®\k')(k\, (3.22) 

jkj'k' 

and the Schmidt number of U is given by the rank of the 
matrix Ujji y kk' = Ufk'jk- But Uj>k',jk = Vj'k'ojka, so 



up to reordering of the columns V = [U\ ■ ■ ■]. It follows 
that the rank of V is at least as great as the rank of U, 
and thus Sch(U) > Sch([7). Taking logarithms of both 
sides we get Kn ar (V) > Ka^U), which is the reduction 
property. 



3. Application to the log-rank lower bound 

As an illustration of the power of the framework we 
have just developed, we now apply it to the analysis of 
a computational problem of considerable interest: the 
communication cost of a distributed computation. 

We consider two separate problems in distributed com- 
putation, the first related to distributed computation of 
a classical function, the second to distributed computa- 
tion of the quantum Fourier transform. The first prob- 
lem may be stated as follows. Suppose Alice and Bob 
are in possession of classical data strings x and y, respec- 
tively. They wish to compute some joint one-bit function 
f(x,y) of the data strings. To accomplish this task they 
are only able to do arbitrary local quantum operations 
and to communicate qubits. This is the key problem of 
quantum communication complexity [45, 46]. 

One of the major results in the fields of quantum 
and classical communication complexity is the log-rank 
lower bound. This states that the minimum number of 
bits (or qubits) of communication required to compute 
f(x,y) is bounded below by log[rank(— l^ x,v ^)], where 
(— l)/( x 'f) is the (x, y)th entry of the communication ma- 
trix. Mehlhorn and Schmidt [47] proved this result for 
classical communication complexity. The log-rank con- 
jecture of communication complexity [48] states that, up 
to a polynomial factor, the log-rank lower bound is satu- 
rated, that is, there is a protocol to compute f(x, y) using 
polylog[rank(— \f( x >y))~\ bits of communication. 

Although quantum protocols are potentially more pow- 
erful than classical, it was pointed out by Buhrman, 
Cleve and Wigderson [49] that [45, 46] contains an im- 
plicit proof of the log-rank lower bound in the quantum 
case. This result was extended to the model in which 
pre-shared entanglement is allowed by Buhrman and de 
Wolf [50]. 

The framework introduced above and the results we 
have proved about the Hartley strength allow us to give 
an almost trivial proof of the log-rank lower bound in the 
case when only qubit communication is allowed, with no 
pre-shared entanglement. The proof is as follows. Sup- 
pose we have a protocol in which Alice and Bob com- 
pute f(x, y) using k qubits of communication. Then it 
is not difficult to see that they can also compute f(x,y) 
using at most k swap gates and no qubit communica- 
tion. Using Bennett's techniques [51] of reversible com- 
putation, the protocol may be modified (using only lo- 
cal unitary operations) to give what Cleve et al. [52] 
called a clean protocol effecting the unitary transforma- 
tion \w A }\x}\y)\w B ) — > (-l) f{x ' v ' ) \w A )\x)\y)\wB), where 
\wa) and \wb) are local work qubits for Alice and Bob. 
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The clean protocol uses only 2k SWAP gates. Let V 
be the unitary effected by the clean protocol, and let 
U\x)\y) = (—l)t( x ' y '\x,y). Then by the reduction prop- 
erty followed by the chaining property we have 

Kh.AU) < K Kai (V) < 2fcX H ar(swAp) - 4fc. (3.23) 

But u = Y. xy (- l ) f(x " v) \ x )( x \ ® \v){y\ from which it; fol - 

lows that Sch(f7) = rank((-l)/( x '«)). Combining this ob- 
servation with Eq. (3.23) gives the log-rank lower bound 



k > j log 



rankf-l^**) 



(3.24) 



The second problem in distributed computation we 
consider is the distributed computation of a unitary oper- 
ation such as the quantum Fourier transform U on m + n 
qubits (m < n), where Alice is in possession of the first 
m qubits, and Bob is in possession of the remaining n 
qubits [76]. How many qubits of communication must 
Alice and Bob do to compute Ul Suppose it is possible to 
achieve it with just k qubits of communication. Then, as 
in the discussion of the log-rank lower bound, it must also 
be possible to implement the quantum Fourier transform 
in a model in which no qubit communication is allowed, 
but in which Alice and Bob can apply k swap gates to 
their qubits. Applying the reduction and chaining prop- 
erties we conclude that Kn ar (U) < fci^Har(sWAP), and 
thus we obtain the lower bound k > 2m, which agrees 
with the m = n result obtained in [4] . 



IV. THE CANONICAL DECOMPOSITION 

Before we describe our results about measures of dy- 
namic strength, we pause to explore a useful representa- 
tion theorem for two-qubit unitary operators, the canon- 
ical decomposition of Khaneja, Brockett, and Glaser [53] 
(see also Kraus and Cirac [26] for a simple, constructive 
proof). This decomposition is an extremely valuable tool 
which characterizes the non-local properties of any two- 
qubit unitary with only three parameters, 9 X , 9 y , and 
9 Z [77]. For appropriate one-qubit unitaries A\, A 2l B\, 
and £>2, 

where — < 9 a < j . For convenience, define the canon- 
ical form of U to be U = (A[ <g> B\)U(A\ <g> b\); up to 
local unitaries, U is equivalent to U. 

Since X <g> X, Y <g> Y, and Z <g> Z all commute, we may 
expand U as 

U = (c x I <g> J + is x X <g> X) x (c y I <g> J + is y Y <g> Y) 

x(c z I®I + is z Z®Z), (4.2) 

where c a = cos(# Q ), s a = sm(8 a ). Multiplying the ex- 
pression out yields 

C X Cy c z + is x s y s z )I <g> I + (c s z + is x c y c z )X (g> X 
+ (s x c y s z + ic x s y c z )Y <g> Y + {s x s y c z + ic x c y s z )Z ® Z. 

(4.3) 



This expression is essentially in Schmidt form: up to a 
constant the Schmidt coefficients are just the magnitudes 
of the coefficients appearing in front of the four terms. 
Eq. (4.3) enables us to deduce the following result [78]: 

Proposition 3 There exist two-qubit unitary operators 
with Schmidt number 1, 2, and 4, but not 3. 

This is a surprising result because it reveals unexpected 
structure in the space of two-qubit unitary operators. It 
is tempting to speculate on the existence of similar struc- 
ture for more general unitary operators. We conjecture 
that, in a d x a" system, there exist unitary operators 
with Schmidt number k if and only if k divides dd' . An 
alternative conjecture, which we believe is less likely, is 
that unitary operators with Schmidt number k exist if 
and only if k and dd' are not coprimc. 
Proof: It is straightforward to see that unitaries with 
Schmidt number 1, 2, and 4 exist, so it only remains 
to show that there exist none with Schmidt number 3. 
Suppose U has Schmidt number 3. Then the canoni- 
cal form of U, U , must have exactly one of the terms 
in Eq. (4.3) equal to zero. Without loss of generality, 
suppose the J ® I term is zero. (If, for example, the 
X <g> X term is zero, then we multiply by X ® X to ob- 
tain a unitary with the I ® I term zero.) Then we must 
have c x c y c z — s x s y s z — 0, and therefore c a = for 
at least one value of a, and sp = for at least one 
value of p. Note that a cannot be equal to /3 since 
s a + c a = 1- By symmetry it suffices to assume that 
a is x and [3 is y, in which case we obtain a unitary of 
the form is x c v c z X® X + s x c y s z Y(g>Y, which has Schmidt 
number at most 2. ■ 
Now suppose that U has Sch(f7) < 2. Then, up to lo- 
cal unitary operations, it has the form of Eq. (4.3), with 
exactly two of the terms non-zero. As mentioned in the 
previous proof, we can always ensure that the I® I term 
is non-zero. Furthermore, conjugating by local unitaries, 
we can ensure that the other non-zero term is X ® X. 
Thus, up to local unitary equivalence, U has the form 
U = aI®I+bX®X, for some non-zero a and b. Further- 
more, we may assume that a is real, since we can multiply 
U by the local unitary operation (e^I) (g> /. Unitarity of 
U then implies that 

I(g>I = U i U= (a 2 + \b\ 2 )I(g>I + a{b*+b)X(g>X, (4.4) 

from which we deduce that a 2 + |&| 2 = 1 and a(b* +b) = 0. 
Since s / 0, i must be pure imaginary. Thus we have 
a = y/l — p, b = iy/p for some < p < 1. We have 
proved the following: 

Proposition 4 Let U be a two-qubit unitary operator 
with Schmidt number 2. Then, up to local unitary equiv- 
alence, U has the form: 



U = y/l-pl® I + i^/pX ® X. 



(4.5) 



11 



V. STRENGTH MEASURES BASED ON 
ENTANGLEMENT GENERATION 

In this and the following section we explore some of the 
strength measures defined by us and other authors, not- 
ing relations between them, and connections to our ear- 
lier operational questions. We also prove several results 
about which measures obey which axioms/properties, 
summarized in Table I at the end of this paper. 

We start in this section with strength measures based 
on entanglement generation. More is known about these 
measures because they use the relatively well-developed 
field of state entanglement. It seems likely to us that, 
although these are natural measures to consider first, in 
the long run they may not be the most useful. Since they 
are based on static resources, they may not provide much 
more insight when applied to dynamics. We consider two 
classes of entanglement-based strength measures: the en- 
tanglement generating capacities of quantum operations 
without initial entanglement, and entanglement generat- 
ing capacities with the possibility of initial entanglement. 



A. Entanglement generation without prior 
entanglement 

Recall the definition of K E {U), Eq. (3.1): K E {U) = 
maxi a \ im E{U\a)\(3)). K E {U) measures the maximum 
amount of entanglement generated by a single application 
of the unitary U without initial entanglement. We show 
that K E and if sch are related to each other in interesting 
ways: (1) i^sch is a lower bound for K E ; and (2) K E is 
equal to i^sch for a class of two-qubit unitaries. We also 
give some numerical evidence demonstrating that K E is 
not equal to i^sch for certain unitaries; see Fig. 2. To 
make this discussion easier, we begin by discussing of the 
properties satisfied by K E and i^sch, including a demon- 
stration of the striking property that K E is superadditive, 
that is U ® U can sometimes generate strictly more than 
twice as much entanglement as U alone. Finally, we ex- 
tend the definition of K E and i^sch to general quantum 
operations, and prove that K E > Kgch still holds. 

1. Properties of Ke and Ks ch 

Beginning with the three axioms, it is easy to see that 
both K E and i^sch satisfy non-negativity, locality, and lo- 
cal unitary invariance. (As we have only defined K E and 
-ft'sch for unitaries, the axioms and properties we discuss 
here are restricted to this case.) 

We now turn to the properties of #sch, which are very 
similar to those of i^Har- i^Sch clearly satisfies the proper- 
ties of exchange symmetry, time-reversal invariance, and 
stability with respect to local ancillas, since none of these 
operations change the Schmidt coefficients. The argu- 
ment that i^sch is continuous is slightly complicated, and 
will be described in the next paragraph, i^sch is strongly 



additive, i.e. K Sch (U <g> V) = K Sch (U) + K Sch (V). To see 
this, recall that if U and V have Schmidt decompositions 
u = Ej s j A j ® B j and V = Efc*fcCfc <g> D k , with Aj, 
Bj, Ck and Dk acting on systems Ai, A2, B\, and £> 2 , 
respectively, then the Schmidt decomposition of U <g> V 
with respect to A1B1 : A2B2 is given by Eq. (3.18): 

U®V = ^2 Sjt k {Aj <g> C\) <g> {Bj <g> D k ). 

jk 

Using properties of the Shannon entropy we find that 

AWl/ 8 F) = tf ({H}) 

-»({&})+»({&}) (51 » 
= K Sch {U)+K Sch (V). 

To see that Kg c h is continuous, expand 

U= E U 3j'M'\i)(k\®\f)(k'l (5-2) 

jj'kk' 

where the comma separates row and column indices. 
Since \j)(k\ and \j'){k'\ are orthonormal operator bases, 
it follows that the Schmidt coefficients of U are 
just the singular values of the matrix U defined by 
Ujkj'k' = Ujji.kw- Consider the matrix norm ||A|| = 
max|^ where the maximization is over unit vec- 

tors \tp). Ksch is a continuous function of the Schmidt 
coefficients, and the Schmidt coefficients are continuous 
functions of the matrix U, with respect to matrix norm. 
This follows from the fact that the singular values of a 
matrix are continuous in the matrix (see, e.g., Chapter 3 
of [54]). Thus ifsch is a continuous function of U with 
respect to the matrix norm. 

We have demonstrated numerically that i^sch does not 
satisfy chaining; see Fig. 1. 



i^sch also violates the reduction property. To see this, 
suppose a Toffoli gate V is applied to three qubits ABC, 
with A acting as the target qubit. Suppose C is initially 
prepared in the |1) state, so V|^)|l) = (£7|^>))|1), where 
U is the CNOT gate, and \ip) is an arbitrary state of AB. 
It is not difficult to verify that K Sch {U) = #(1/2, 1/2), 
while K Sch {V) = #(1/4,3/4), so K Sch {V) < K Sch (U), in 
violation of the reduction property. 

The properties of K E are somewhat more difficult to 
elicit. K E is easily seen to satisfy the exchange sym- 
metry property. Numerical studies of the time-reversal 
invariance property have been inconclusive, although we 
speculate that for two-qutrit unitaries time-reversal in- 
variance will not be obeyed. The discussion of conti- 
nuity is somewhat complicated, and is described in the 
following paragraph. K E is stable with respect to ancil- 
las, since it already allows for the possibility of arbitrary 
ancillas. It is also easy to see from the definition that K E 
satisfies the reduction property in the same sense that 
the Hartley strength satisfies the reduction property. 
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FIG. 1: Numerical violation of the chaining property for Ks c h- 
U and V are two-qubit unitaries chosen by first generat- 
ing random unitaries, and then using a Nelder-Mead sim- 
plex minimization algorithm to prepend and append local 
unitaries to generate U and V maximizing the violation of 
Ks c h(UV) < Ksch(U) + Ksch(V). If ifsch satisfied chaining, 
then all the points (corresponding to pairs of unitaries U and 
V) would lie on or below the line. 

We now outline a proof that K E is continuous. To 
prove this, we need to introduce a metric on the space of 
unitary matrices. We use the matrix norm, D(U, V) = 
\\U — V\\ = max|^ ||(I7 — V)|^)||, where the maximum 
is over all unit vectors \ip). Choose \a), \[3) such that 
K E (V) — E(V\a)\/3)). Our earlier discussion shows that, 
without loss of generality, we may assume \a) lives in a 
(i^-dimensional space, and \(3) lives in a rig-dimensional 
space. It follows from the definition that 

K E (U) > E(U\a)\f3)) (5.3) 

The results of [42] (see also [55]) imply that, provided 
|| |V>- |0>II < 1/6, 

\E(\i > ))-E(\ct>))\<2\\\i > )-\<t>)\\log(dcdv) 

+ 7/(211 |V> - |0)||), (5.4) 
where r](x) = —xlog(x). Thus, provided \\U — V\\ < 1/6, 

\E(U\a)\(3)) E{V\a)\[3))\ < 2\\U - V\\ log(d A d 2 B ) 

+ V (2\\U-V\\). (5.5) 

Combining this result with Eq. (5.3) and the fact that 
K E (V) = E(V\a)\(3)), we obtain 

K E (U) > K E (V) -2\\U- V\\ log(44) - V (2\\U - V\\). 

(5.6) 

By symmetry the same inequality holds with U and V 
interchanged, and thus 

\K E (U)-K E (V)\<A\\U-V\\log(d A d B )+ V (2\\U-V\\) 

(5.7) 



whenever \\U — V\\ < 1/6, which is the desired continuity 
equation. 

What about the additivity properties of K E 1 Intu- 
itively, we expect the amount of entanglement generated 
by two copies of U is no greater than twice the maximum 
generated by one use of U. However, this intuition fails 
when ancillas are allowed. We show below that, unlike 
^Sch, K E is superadditive. The proof requires some facts 
about the relationship between K E and A'sch, so we prove 
this result at the end of Sec. V A 2. Since K E is stable 
with respect to local ancillas, subadditivity of K E and 
Proposition 2 imply that K E does not satisfy chaining. 

2. Relations between Ke and Tfsch 

In this subsection, we explore some relations between 
K Sch and K E . 

Lemma 5 For all unitaries U, Ks c h(U) = E(U\a)\(3)) 
where \a) is a maximally entangled state of system A 
with an ancilla TZa, and \0) is a maximally entangled 
state of system B with an ancilla 1Z& . 

Proof: Let A and B label Alice's and Bob's systems, re- 
spectively. Alice introduces an ancilla TZa that is a copy 
of her system. She prepares A and TZ A in a maximally 
entangled state, \a) = where d A is the 

dimension of system A (and hence also of system 1Zj[)- 
Bob does the same thing, preparing \(3) = J+ |j)|j), 
where <ig is similarly the dimension of B. 

Let U — J2 t siAi®Bi be the Schmidt decomposition of 
U (Eq. (2.1)). Alice and Bob apply U to AB, obtaining 

U\a)\(3) = s^B^) = -^=\ai)\h), (5.8) 

where we define \a{) = ^fdA~Ai\a) and \bi) = \fd^ Bi \ f3) . 
\ai) and \bi) are orthonormal bases. For example: 

(a k \ai) = d A (a\A\Ai\a) = tr^Aj = 8 U . (5.9) 

Therefore, U\a)\0) has entanglement H ({sf / (d^dg)}) 
which is equal to Ks c h{U). ■ 
From this lemma, it follows that K E (U) is bounded 
below by Ks c h(U). We also show that they are equal for 
certain two-qubit unitaries: 

Theorem 6 K E (U) > Ks c h(U) for all unitaries U. 

Theorem 7 K E (U) = K~s c h{U) for all two-qubit uni- 
taries U satisfying Sch(J7) < 2. 

Proof: When Sch(J7) = 1, U is a local unitary and hence 
K E (U) = K Sch {U) = 0. 

Suppose Sch(J7) = 2, in which case Proposition 4 im- 
plies that U may be expanded as 

U = {A l ® B x ) (y/l-pl® I + iy/pX O a) (A 2 ® B 2 ). 

(5.10) 
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Let U = \/l — p I® / + i^fpX ® X. We have seen in the 
previous section that Ke and ifsch are both invariant 
under local unitaries, so we have Ke{U) — Ke(U) and 
K S ch(U) = K Sch (U). 

We can calculate Ks c h(U) and Ke{U) directly. 
Ksdi{U) is equal to H(l—p,p) = H(p), the binary Shan- 
non entropy. To calculate Ke{U), we substitute U into 
the expression Eq. (3.1) for Ke, giving 



K E (U) = max S 



[l-p)\a)(a\+pX\a)(a\X 



+ t^/p(l-p)(/3\X\(3)(X\a)(a\ ~ \a)(a\X) 



(5.11) 



where S is the von Neumann entropy, and its argument 
is a state of ATZ A . Now we use the fact that a projec- 
tive measurement on A1Z A cannot decrease its entropy 
(see Chapter 11 of [5]). We measure in an orthonormal 
basis containing the elements \a) and |a_i_), where |aj_) 
is chosen so that, up to an unimportant global phase, 
X\a) — cos (f>\ a) + sin</>|aj_) for some 4>. We obtain 

K E (U) < max5 [(1 - p)\a){a\ + p(a\X\a) 2 \a)(a\ 

\a) 



+ p\(a ± \X\a)\ 2 \a 1 _)(a ± \] 
< maxiJ(l — p + p cos 2 <p,psm 2 (f>). 



(5.12) 



If p < i, the maximum occurs for <f) = § and Ke{U) < 
if (p) = K Sch (U). (If p > i, apply X (g> X to U to swap 
the role of p and 1— _p.) Since, by Th. 6, K E (U) is greater 
than or equal to Ks c h(U), we must have equality. ■ 
We show below that Ke is superadditive while -ftTsch is 
additive, which implies that they are not equal for cer- 
tain unitaries. We have also shown this numerically by 
calculating both functions for a particular class of uni- 
taries, the Schmidt number 4 family parametrized by p, 
denoted U p in Eq. (2.5). Fig. 2 plots both K E (U P ) and 
Ks c h(U p ) as a function of p, and also their difference. 



We now have the tools required to prove that Ke is 
superadditive, as promised at the end of the last section. 



Theorem 8 Ke is 

taries U such that 



superadditive, i.e. there exist uni- 



K 



A 1 A 2 .B 1 B 2 



AiBi 



U A2B2 )>2K^(U AlBl ). (5.13) 



where the subscripts on U indicate the subsystems to 
which it is applied. 



Proof: Let U = ^l-pl ® I + iy/pX <E> X. We show 
that additivity is violated for certain values of p. (We 
will only add subscripts where necessary.) 

Since U has two Schmidt coefficients, Th. 7 implies that 
Ke(U) = Ksch(U). Therefore, the right-hand side of 
Eq. (5.13) is 2K E {U) = 2K Sch {U) = 2H(p). 

To obtain the violation of additivity Eq. (5.13) we now 
construct specific states |a) and \ j3) of A and B for which 
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FIG. 2: Plots of K E {U P ) (dots) and K Sc h(U p ) (solid) as func- 
tions of p, and of the difference Ke(U p ) — Ks ch (U p ) (dashed), 
demonstrating that Ke(U p ) 7^ Ksch(U p ) for some values of p. 



E(U\a)\P)) > 2H(p). To do this, we apply £7® U to two 
pairs of systems, as depicted in Fig. 3, where we have 
omitted the ancillas as they turn out not to be necessary 
for our construction of \a) and \(3) . Let |a) = ( 1 00) + 
\ll))/^/2 be a state of Alice's system A1A2 and \(5) = 
(|00) + |H))/V2 be a state of Bob's system BiB 2 - 




FIG. 3: Diagram of U <giU applied to systems A1B1 and A-2.B2- 
Note that A1A2 starts out in the maximally entangled state 
I a), and B1S2 starts out in the maximally entangled state \/3), 
so A1A2 is not initially entangled with B\B2- 



We make use of a handy identity to calculate 
E(U\a}\(3}). Since |a) and \(3) are maximally entangled, 
a calculation shows that for any two-qubit unitary U , 



(U 



AtBt 



I A 2 B 2 ) I ") I/?) = ( I A 1 



Si 



Ul 2 B 2 



(5.14) 
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where the transpose is taken in the basis 
{|00), |01), 1 10), |11)}. This is illustrated in Fig. 4. 




FIG. 4: Illustration of the identity Ua 1 b 1 ® lA 2 B 2 \a:)\f3} = 
lA lBl ®Ul 2B2 \a)\l3). 



For the unitary we are considering, U T = U, so that 
Eq. (5.14) implies 

E(U AlBl ® U A2B2 \a)\(3)) = E(I AlBl ® U 2 MB2 \a)\(3)). 

(5.15) 

We may now apply Lemma 5, considering A\ and B\ 
as the ancillas to Ai and B2, respectively. We see that 
E(U AlBi m A2 B 2 \a)\f3}) = K Sch (U 2 ). Observing that U 2 
is a unitary with two Schmidt coefficients, 

U 2 = (1 - 2p)I ®I + 2iy/p(l -p)X®X, (5.16) 
we obtain 

K E (U ®U)> E(U ® U\a)\[3)) = H[(1- 2p) 2 ] , (5.17) 

so we have reduced the problem to showing that there 
exist values of p such that H [(1 - 2p) 2 ] > 2H(p). The 
existence of such values is shown in Fig. 5 [79]. ■ 



3. Extension to general quantum operations 

Our results to this point have primarily concerned 
strength measures for unitary operations. In this sub- 
section, we obtain some results for general quantum op- 
erations, proving generalizations of Lemma 5 and Th. 6 to 
quantum operations. We will not do a detailed investiga- 
tion of the axioms and properties satisfied by these mea- 
sures for general operations, although arguments similar 
to the unitary case mostly go through. 

The first step is to generalize our definitions of Ke 
and ifsch- In order to generalize K E (Eq. (3.1)) to quan- 
tum operations, we must choose an entanglement mea- 
sure which applies to mixed states as well as pure states. 
We use the entanglement of formation [56]: 

F(p) = min^.^(|^)) (5.18) 




0.2 0.4 0.6 0.8 1 

P 



FIG. 5: A plot of 2H(p), H [(1 - 2pf] , and their difference. 

where the minimization is over all pure state decom- 
positions {pj, \ipj)} of P, and E is the entanglement 
of pure states. Note that any two decompositions of 
p are related by a right unitary matrix Ujk- p = 
Z)jPj#j><^jl = 2~2k1k\(f>k)((f>k\ if and only if [57-59] 
y/Pjl^j) — J2j Ujky/qk\<f>k)- We take as our generalized 
Ke{£) the maximum entanglement generated by £ over 
all separable input states pa-.b- 

K E (£)=ma X F[£(p A:B )}. (5.19) 

PA-.B 

Note that F o £ is a convex function maximized on the 
convex set of separable states, {pa-.b}, an d therefore 
F [£(pa-.b)} achieves its maximum for extreme points of 
the set of separable states, i.e. pure product states. 

To generalize -Ksch, let £ be a quantum operation with 
operation elements {Gfe}: £{p) = J2k GkpG\. £ can be 
decomposed differently as £{■) — ^ . Fj ■ Fj if and only 
if [60] the two sets of operation elements are related by 
a right unitary matrix: Fj = ^ fc UjkGk- By analogy 
with the entanglement of formation, a natural definition 
of K Sch (£) is 

K Sch (£) = min V ^f^ Ks^Fj), (5.20) 

3 

where Ksch(Fj) is given by Eq. (2.15), and the minimiza- 
tion is over all possible decompositions of £ into opera- 
tion elements. The coefficients ti(F^ Fj) / '(d A d B ) form a 
probability distribution. A physical interpretation is as 
follows: if Ksch(Fj) is the strength of the operation Fj, 
then Ksch(£) is the expected strength of £, minimized 
over all possible decompositions of £ . 

First, we prove two lemmas generalizing Lemma 5. For 
the remainder of this section, let |a) be a maximally en- 
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tangled state of system A with an ancilla TZ_a, and \(3) be 
a maximally entangled state of B with an ancilla TZ^. 

Lemma 9 For all operators Q, 

Ks*{Q) = E (y/$foQ\ a )\0)) . (5.21) 

Proof: Recall that K Sch (Q) = H ({sf /tr(Q+Q)}) , so 
we need only calculate the right-hand side of Eq. (5.21). 
Expand the state Q\a}\[3} as 

Q\a)\0) =Y J *iMa)Bi\P) = £ -^=h)h), 

(5.22) 

where Y^i s iAi <S> Bi is the Schmidt decomposition for Q, 
and \ai) = \/d^ Ai\a) , \bi) = \fds Bi\f3) are orthonormal 
bases for their respective systems. The result follows. I 

Lemma 10 For any quantum operation £, let a = 
£{\a)(a\ <g> \/3)(f3\). Then K Sch (£) = F(a), where F is 
the entanglement of formation. 

Proof: Let Fj be the set of operation elements for £ 
achieving the minimum in the definition of i^sch- Then, 
applying the definition and Lemma 9, we have 

(5.23) 

Noting that 

is an ensemble for cr, we deduce that Ks c h(£) > F(a). To 
prove the reverse inequality, suppose o = ^2 k Pk\4>k){<j>k\ 
is the minimizing decomposition for the entanglement of 
formation of cr. Note that a can also be decomposed as 

<j = J2 j F 1 (\u)(a\®\[3}(f3\)Fj. (5.25) 

The minimizing decomposition is related to the decom- 
position from Eq. (5.25) by a right unitary matrix U: 
\fPk\ ( t ) k) = J2jUkjFj\a}\[3). This unitary freedom is 
identical to the freedom in the operator-sum decompo- 
sition, so the set of elements Gk = X)j Ukj Fj is also 
an operator-sum decomposition for £ , as well as giving 
the minimizing decomposition of cr, that is y/p~k\ <j>k) — 
Gk\a)\(3). This gives us the desired inequality, 

> K Sch (£). 



The desired bound on Ke now follows: 

Theorem 11 Ke{£) > K^ c h{£ ) for all quantum opera- 
tions £. 

Proof: The result follows immediately from the previous 
lemma and the fact that 

K E {£) = m & xF[£( PA .. B )} >F[£(\a)(a\ ® . 

(5.27) 



B. Entanglement generation with prior 
entanglement 

In this section we consider the largest change in en- 
tanglement which can be caused by a unitary U , us- 
ing both ancillas and prior entanglement, as defined in 
Eq. (3.2) and repeated here for convenience: K\e{U) = 
sup^ \E(U\ip}) — E(\ip))\, where U acts on the combined 
system AB, and \ip) is an arbitrary state of AB plus their 
ancillas, 1Z A and 7?.g. We show that, although K AE in- 
volves a more difficult maximization than Ke, and may 
therefore be more difficult to work with, it satisfies more 
of the axioms and properties described in Sec. Ill B than 
Ke does. Incidentally, since K A e and Ke have different 
properties they can not, in general, be equal. 

We first show that K ae obeys the three axioms. K A e 
is clearly non-negative and satisfies local unitary invari- 
ance. To show that K\e satisfies locality is only slightly 
more involved. If U = A® B, then K AE (A <g> B) = 
sup w \E(A <g> B\i[))) - E(\ip))\ = 0. On the other hand, 
since K A e(U) > Ke(U) and we know that Ke{U) sat- 
isfies locality, K\e{U) — only if Ke(U), which implies 
that U is a local unitary, as required. 

Second, we show that K A e satisfies Properties 1, 2, 
4, and 6-8. Properties 1 and 2, exchange symmetry and 
time-reversal invariance, are easily seen to be true. We 
do not know whether property 3, continuity, is satisfied. 
The argument used to establish that Ke is continuous 
does not work in this instance, because we do not have 
any bound on the size of the ancilla that A and B may 
use. If such a bound could be established then a similar 
continuity bound to that used for Ke could be proved. 
Next, we show that K A e obeys chaining, Property 4. For 
any two unitaries U and V, 

K AE (UV) - sup \E(UVm - £(|V»| 

IV-) 

= sup \E(UVm - E(V\t/>)) + E(vm - 
IV-) 

< sup \E(U\4>))-E(\4>))\+sup\E(V\iP))-E(\iP))\ 
\4>)=VW |V) 

= K AE (U) + K AE (V). (5.28) 

Property 6, stability with respect to ancillas, holds since 
K A e already allows the possibility of arbitrary ancillas. 
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Therefore, by Proposition 2, Kae also satisfies strong 
subadditivity, Property 8. Finally, we note that the defi- 
nitions immediately imply that K^e satisfies the reduc- 
tion property Property 9. 

VI. STRENGTH MEASURES BASED ON 
METRICS 

In this section we consider a class of strength measures 
motivated by the axiomatic approach. This is in con- 
trast to Sec. V, where we studied strength measures based 
on entanglement generation. The strength measures we 
study here are based on metrics. We explore the axioms 
and properties obeyed by these measures when different 
constraints are placed on the underlying metrics. We 
derive an exact, analytic formula for one particular mea- 
sure. Finally, we examine the potential of these measures 
for analyzing quantum computational complexity, as de- 
scribed in Sec. Ill B. 

Recall the definition of a metric. Let X be a set. A 
metric is a real function D:IxX^I satisfying the 
following properties for any three elements x, y, z of X: 

(1) D(x, y) > with equality if and only if x = y 

(2) D(x,y) = D(y,x) (symmetry) 

(3) D(x, z) < D(x, y) + D(y, z) (triangle inequality) 

Given a metric D, the corresponding strength measure 
Kd{U) is the minimum distance between U and the set 
of local unitaries LU: 

K D (U) = min D(U,L). (6.1) 

LeLU 

The set LU varies depending on context. The most com- 
mon case is where U is a two-qudit unitary acting on the 
space AB and LU is the set of products of one-qudit uni- 
taries, Krj(U) = mmA,B D(U, A® B). Analogues of the 
definition of Kr> were introduced to quantify state en- 
tanglement by Vedral et al. [37], and have been studied 
in considerable detail, proving to be a fruitful approach 
to quantifying state entanglement. 

More generally, if U acts on a composite of systems, 
Ai,A2,---,A m , there are several notions of "local", 
which we differentiate with superscripts. For example, 
suppose U acts on ABC. One notion of "local unitary" 
corresponds to unitaries of the form A ® B ® C, so that 
Kg B:C (U) = min.A,B,c D(U, A ® B ® C). A different 
division into subsystems leads to a different measure: 
Kfi BC (U) = min4,B D(U, A ® B), where A acts on sys- 
tem A but now B is any unitary on BC. 

A. Properties of strength measures based on 
metrics 

One reason for studying strength measures based on 
metrics is that the properties of the strength measure 



may be controlled by varying the properties of the un- 
derlying metric. We consider strength measures based 
on: (1) arbitrary metrics; (2) metrics invariant under lo- 
cal unitaries; and (3) metrics invariant under any unitary. 
Each extra requirement causes the strength measure to 
obey extra axioms and properties from Sec. Ill B. Since 
we know of no general way to characterize families of 
metrics, in this section we do not consider any of the 
properties applying to families (Properties 5-9). There- 
fore, throughout this section we assume Kb = K^' B . 

The metric properties are easily seen to guarantee that 
the axioms of non- negativity and locality hold for all Kn . 
An elegant fact is that the metric properties alone also 
imply that Kb satisfies the continuity property: 

Lemma 12 For any two unitaries U and V, and any 
metric D, \K D (U) - K D (V)\ < D(U, V). 

Proof: Choose A and B such that K D {V) = D(V, A ® 
B). By definition K D (U) < D{U,A®B) 1 and by the tri- 
angle inequality D(U, A® B) < D(U,V) + D(V, A® B) = 
D(U, V) + K D (V). Thus K D (U) < D(U, V) + K D (V), 
which may be rearranged to give Ke>{U) — Kd(V) < 
D(U, V). By symmetry, K D {V) - K D {U) < D(U, V). ■ 

If D is locally unitarily invariant, i.e. , D(U, V) = 
D[(A®B)U,(A®B)V] = D[U(A® B),V(A® B)}, then 
A'psatisfies local unitary invariancc. 

Finally, suppose the metric satisfies full unitary invari- 
ance, so that D(U, V) = D(WU, WV) = D(UW, VW) 
for any unitary W. Then Kd satisfies two additional 
properties. The first is exchange symmetry, which is 
easily proved. The second is chaining, Kr>{UV) < 
K D (U) + K D (V). To see this, suppose A® B and C® D 
minimize Kd(U) and Kn(V), respectively. Then 

K D {UV) < D [UV, {A ® B){C ® £>)] 
< D [UV, U{C ®D)]+D [U{C ®D),{A®B){C® D)] 
= D(V, C®D) + D{U, A® B) 

= K D (U) + K D (V). (6.2) 

B. An explicit formula for the Hilbert-Schmidt 
strength of a two-qubit unitary 

In this section we consider an example of a metric- 
based strength measure, the Hilbert-Schmidt strength 
A"hs induced by the unitarily invariant Hilbert-Schmidt 
norm on operators, ||Q||hs = \/tr(QtQ). More explicitly, 
for a bipartite unitary operation U we define 

A-Hs(^) = min||[/-yl®B||HS, (6.3) 

A..B 

where A and B are local unitary operators on the respec- 
tive subsystems. We now exhibit an explicit formula for 
the Hilbert-Schmidt strength in the two-qubit case. 

The statement of the result is simplified by first making 
some definitions and observations. Let U be a two-qubit 
unitary operation with canonical decomposition 

U = {A l ®B l )e^ x ® x+e - Y ® Y+ ^ z ® z \A 2 ®B 2 ). (6.4) 
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Because of local unitary invariance the Hilbert-Schmidt 
strength depends only on the parameters 6j, that is, we 
can ignore the local unitary operations A\p and B\^- 
Without loss of generality we assume U is in canonical 
form, that is, A\ = B\ = A 2 = B 2 = I- 

We define |0 O ) = (100) + |H))/\/2, and \<j>j) = (I ® 
<jj)\<po) for j = 1,2,3 where we write ctq , o~\ , o 2 , 03 to 
denote I, X, Y, Z. Note that the set \<j) 3 ) for j = 0, 1, 2, 3 
is the Bell basis, up to phases. A simple but tedious 
calculation verifies the useful formula (<pj\a k ® o~i\4>j) — 
SkiHjk, where the 4x4 matrix H is 



H = 



The H matrix can also be used to evaluate the eigen- 
values of U. Because X ® X, Y ® Y, and Z ® Z are 
diagonal in the \<pj) basis, U may be written in diagonal 
form as U = J2j ^j\ ( f ) j)( ( l ) j\> where Aj are the eigenvalues 
of U. These eigenvalues are evaluated as follows: 



,i(8iX®X+0 2 Y®Y+0 3 Z®Z) 



= exp ^i'y]0k((f>j\o'k ® °~k\<f>j)J , 



(6.5) 



where in the last line we used the fact that all three a k 1 
(Tfc are diagonal in the \cf)j) basis. Substituting ((fijWk 1 
(Ti\4>j} = 5kiH jk we obtain: 



Xj = exp i^HjkOk 

\ k=l 



(6.6) 



Theorem 13 For a two-qubit unitary U with canonical 
decomposition Eq. (6.4), the Hilbert-Schmidt strength is 
given by the formula 



Kns(U) = 4/8-2 max V . XjH 



0<k<3 



ijk 



(6.7) 



The minimizing local unitary is A®B — e t9 ak®cr k where 
k achieves the maximum in the expression above, and 6 
is the argument o/Ej XjHj k . 

Proof: Simple algebra shows that 

K HS (U) 2 = min[8-23f?(tr[[/ t (A®B)])] , (6.8) 

where Sft(-) denotes the real part. We expand A and 
B in terms of the Pauli operators as A — Efc=o a fc (J fc> 
B = Ef=o (Note that the unitarity of A and B im- 
plies that J2k \ a k\ 2 = E; \bi\ 2 = 1-) Substituting these 
expressions for A and B, and U = J2j gives 



^hs(^) 2 = min 

ak,bi 



8-23? y^X*a k bi{(j)A(T k ® ai\<j> 3 



jki 



(6.9) 



where the minimization is over all a k , bi such that the cor- 
responding A and B are unitary. But (0j |cr fc ® cri\<pj)) — 
SkiHjk, as noted earlier, so this expression simplifies to 

K m {U) 2 = 8 - 2 max 3? (V - fc A*a fc 6 fe ^ fe ) . (6.10) 
The Cauchy-Schwarz inequality implies J2k \ a k°k\ < 1, 



so: 



< max > A*i/ofe 
fe 1^—0' J I 

= max I XjHjk ■ (6-11) 

Equality occurs when a/ = and bi = 5 k ie l6 , where 
k maximizes the right-hand side of the inequality, and 



ei6 Ej X j H . 



ijk 



Ej ^3 H 3k 



. This corresponds to A 



B = e l0 a k ® cjfe, and as described in the statement of 
the theorem. ■ 



C. Applications to computational complexity 

We have seen that strength measures based on unitar- 
ily invariant metrics satisfy many desirable axioms and 
properties. It is natural to ask whether these measures 
might be useful in answering questions about computa- 
tional complexity, as described in Sec. IIIB. 

In order for a family of measures {Kp} to be useful in 
this context, we require {Kr>} to be stable under addition 
of systems (for the remainder of this section, we simply 
write "stable" for this property). This is to ensure that 
the strength of a CNOT gate is independent of the number 
of qubits in the system being studied. It is tempting to 
consider a family of measures {Kr>} whose underlying 
family of metrics is stable, in the sense that D(U, V) = 
D(U ®I,V®I) for any unitaries U and V. However, we 
show here that such metrics give rise to trivial bounds on 
computational complexity. Denote by U a unitary acting 
on n qubits, and let and / be the zero and identity 
operator, respectively, on n qubits. For any such unitary, 



K D (U) 



= min D(U, A 1 ® . . . <g> A n ) 
< ^min [D(U,0) + D(0,A 1 ® 



= 2£>(J,0), 



An)} 

(6.12) 



where to obtain the last line we used the unitary invari- 



ance of D. But I = h® I 2 ® 



I n , where Ij is the 



identity on the jth qubit, so by the metric stability prop- 
erty K D (U) is always bounded by 2£>(J,0) = 2£>(Ji,0), 
which is a constant. Therefore, the lower bound on the 
number of two-qubit gates required to implement a n- 
qubit gate, k > K(U)/K max (Eq. (3.7)), is a constant. 

This shows that any family of metrics which is both 
unitarily invariant and stable cannot give interesting 
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Measure: 


i?Har 


Ksch 


K E 


Kae 


K D itcpLU] K D [U] 


Al Non-negativity 


yes 


yes 


yes 


yes 


yes 


yes 


yes 


A2 Locality 


yes 


yes 


yes 


yes 


yes 


yes 


yes 


A3 LU invariance 


yes 


yes 


yes 


yes 


? 


yes 


yes 


PI Exchange 


yes 


yes 


yes 


yes 


? 


? 


yes 


P2 Time-reversal 


yes 


yes 


? 


yes 


? 


? 


? 


P3 Continuity 


no 


yes 


yes 


? 


yes 


yes 


yes 


P4 Chaining 


yes 


no 


no 


yes 


? 


? 


yes 


P5 System stability 


- 


- 


- 


- 


? 


? 


? 


P6 Ancilla stability 


yes 


yes 


yes 


yes 


? 


? 


? 


P7 Weak additivity 


yes 


yes 


no 


yes 


? 


? 


? 


P8 Strong additivity 


yes 


yes 


no 


yes 


? 


? 


? 


P9 Reduction 


yes 


no 


yes 


yes 


? 


? 


? 



TABLE I: Summary of axioms and properties of strength 
measures, "yes" /"no" indicates whether the strength mea- 
sure obeys the axiom/property. "-" means the property is 
not applicable, and "?" means we do not know whether the 
strength measure obeys the axiom/property. JCd[LU] refers 
to strength measures induced by locally unitarily invariant 
metrics, and .KdRJ] refers to strength measures induced by 
unitarily invariant metrics. 



lower bounds on computational complexity. As noted 
above, unitary invariance is a useful property. On the 
other hand, stability of the family of metrics may not be 
necessary for stability of the induced family of measures. 
So, it may be possible to find a family of unitarily in- 
variant metrics which is not stable, but which induces a 



stable family of measures, and could therefore give useful 
lower bounds on computational complexity. 



VII. SUMMARY AND FUTURE DIRECTIONS 

We have developed the beginnings of a quantitative 
theory of quantum dynamical operations as a physical 
resource. While promising preliminary results have 
been obtained, an enormous amount of work remains 
to be done. (Table I summarizes the properties of the 
strength measures we investigated.) We believe the 
development of this theory offers a concrete path to 
address the fundamental question of quantum com- 
putational complexity: how many one- and two-qubit 
quantum operations are required to do some desired 
quantum operation £1 This will, in turn, allow us to 
answer questions about the relationship of quantum and 
classical complexity classes, and may enable the resolu- 
tion of some longstanding questions in complexity theory. 
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